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To detect changes in two multi-temporal SAR (Synthetic Aperture Radar) more efficiently a new technique is 
discussed in this paper. An unsupervised distribution-free change detection approach is for synthetic aperture radar (SAR) 
images based on an image fusion strategy and a novel fuzzy clustering logarithm. The image fusion technique will be used 
to generate a difference image by using a mean-ratio image and a log-ratio image. DWT will be used for image fusion. 
Wavelet fusion rules based onan average operator and minimum local area energy are chosen to fuse the wavelet 
coefficients for a low-frequency band and a high-frequency band respectively are used in image fusion. A reformulated 
fuzzy local-information C- means clustering algorithm is proposed for classifying changed and unchanged regions in the 
fused difference image. It will enhance the changed information and of reduces the effect of speckle noise. 
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Synthetic aperture radar operates in any weather conditions and provides high resolution images. So it has various 
applications like Surveillance, and Targeting, Foliage and Ground Penetration etc. Image change detection [9] is a process 
that analyses images of the same scene taken at different times in order to identify changes that may have occurred 
between the considered acquisition dates [1]. SAR images suffer from the presence of the speckle noise which degrades the 
quality of image and makes change detection a difficult task. However, SAR sensors are independent of atmospheric and 
sunlight conditions, which make the change detection in SAR images still attractive [1],[5]. 

An unsupervised change detection method [9] is used so we can avoid the requirement of a ground truth image. 
The method comprises of following steps: 1) image pre-processing; 2) producing difference image between the images; 
and 3) analysis of the difference image. The tasks of the first step mainly include co registration, geometric corrections, 
and noise reduction. In the second step, two images are compared pixel by pixel to generate the difference image. For the 
remote sensing images, differencing (subtraction operator) and rationing (ratio operator) are well-known techniques for 
producing a difference image. In differencing technique, changes are measured by subtracting the intensity values pixel by 
pixel between the considered couple of temporal images. 

And in rationing, changes are obtained by applying apixel-by-pixel ratio operator to the considered couple of 
temporal images. However, in the case of SAR images, instead of differencing the ratio operator is typically used operator 
because the image differencing technique is not adapted to the statistics of SAR images and non robust to calibration 
errors [4]. In addition, because of the multiplicative nature of speckles, the ratio image is usually expressed in a logarithmic 
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or a mean scale [2] .With the log-ratio operator, the multiplicative speckle noise can be transformed in an additive noise 
component. Furthermore, the range of variation of the ratio image will be compressed and there by enhances the 
low-intensity pixels while weakening the pixels in the areas of high intensity; therefore, the distribution of two classes 
(changed and unchanged) could be made more symmetrical. However, the information of changed regions that is obtained 
by the log-ratio image may not be able to reflect the real changed trends in the maximum extent because of the weakening 
in the areas of high-intensity pixels. As for the ratio mean operator, the background (unchanged regions) of mean-ratio 
image is quite rough, for the ratio technique may emphasize the differences in the low intensities of the temporal images. 

In the third step, changes are usually detected by applying a decision threshold to the his to gram of the difference 
image. Several thres holding methods [10] have been proposed in order to determine the threshold in an unsupervised 
manner, such as Otsu, the Kittler and Illingworth minimum-error thres holding algorithm (K&I), and the expectation 
maximization (EM) algorithm [6]. It is apparent that this kind of methods requires an accurate estimation of the decision 
threshold. Various fuzzy c-means clustering algorithms [7], [8] were used in change detection like FCM_S, fast 
generalized FCM, FGFCM, and FLICM but all those were not much efficient. FCM was very sensitive to the speckle 
noise, in FCM_S the labelling of a pixel to be influenced by the labels in its immediate neighbourhood. 

However, compared with the original FCM, the computational complexity of FCM_S is significantly increased 
since it computes the neighbourhood term in each iteration step FGFCM significantly reduced the execution time by 
clustering on gray-level histogram rather than on pixels; meanwhile, It is less sensitive to noise to some extent because of 
the introduction of local spatial information. An artificial parameter is applied in their objective functions in order to 
balance between robustness to noise and effectiveness of preserving the details of the image. The selection of parameter is 
not easy to implement since there is no prior knowledge about the speckle noise level. FLICM introduced a term fuzzy 
factor which uses gray-level difference and spatial distance and use of spatial distance is not appropriate in some cases. 
Therefore a new algorithm is introduced. 

In general, it appears clearly from the literature that the whole performance of SAR -image change detection is 
mainly relied on the quality of the difference image and the accuracy of the classification method. In order to address the 
two issues, in this paper, we propose an unsupervised distribution-free SAR-image change detection approach. It is unique 
in the following two aspects: 1) producing difference images by fusing ame an-ratio image and a log-ratio image and 2) 
improving the fuzzy local-information c-means (FLICM) clustering algorithm, which is insensitive to noise, to identify the 
change areas in the difference image, without any distribution assumption. 

This paper is organized into five sections. 

2. PROPOSED METHOD 

Let us consider the two co registered intensity SAR agesXi={xi(i,j),l<i<H,l<j<W}andX2={x 2 (i,j),l<i<H,l<j<W} 
of size, i.e., acquired over the same geographical area at two different times and, respectively. Our objective is aiming at 
producing a difference image that represents the change information between the two times; then, a binary classification is 
applied to produce a binary image corresponding to the two classes: change and unchanged. As shown in Figure 1, the 
proposed unsupervised distribution-free change detection approach is made up of two main phases: 1) generate the 
difference image using the wavelet fusion based on the mean-ratio image and the log-ratio image; and2) automatic analysis 
of the fused image by using an improved fuzzy clustering algorithm. 
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Generate the Difference Image Using Image Fusion 



A discrete wavelet transform (DWT) is used for image fusion. The DWT isolates frequencies in both time and 
space, allowing detail information to be easily extracted from images. Figure 2 shows implementation of a image fusion 
scheme using DWT. First the mean-ratio operator and the log-ratio operator are applied to the two source images used for 
fusion respectively, which are commonly given by 



X m = 1-min (ui/u 2 , U2/U1) 
X 1= llogCX/XOI 



(1) 

(2) 
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Figure 1: Flowchart of the Proposed Change Detection Approach 

Where Uj and |j 2 represent the local mean values of multi temporal SAR images X, and X 2 , respectively. In the 
image fusion scheme based on the wavelet transform wavelet coefficient map is obtained by taking DWT of log ratio and 
mean ratio images. After 
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Figure 2: Image Fusion Technique Based on DWT 

That fusion rules are applied. Here, two main fusion rules are applied: the rule of selecting the average value of 
corresponding coefficients for the low-frequency band, and the rule of selecting the minimum local area energy coefficient 
for the high-frequency band. The fusion rules can be described as follows: 



D F LL =(D m LL +D' L L)/2 



D F f = D m E (i,j),E ,I, E (i,j)<E 1 E (ij)or 



DWj), E' f (i,j) >=E m £ (i,j) 



(3) 



(4) 



Where m and 1 represent the mean-ratio image and the log-ratio image, respectively. F denotes the new fused 
image. D LL stands for low-frequency coefficients. D' E (i,j) represents three high-frequency coefficients at point (i,j) in the 
corresponding sub-images. The local area energy coefficient E E (i,j) can be computed as follows: 



E E (i,j)=2 k£NlJ ID E (k)l 
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Where E E (i,j) represents the local area energy of the wavelet coefficient at point (i,j)in the corresponding 
sub-image, and Ny represents the local window centred on (i,j). D E (k) denotes the value of the k th wavelet coefficient that 
is around the local window. After that inverse DWT is taken to get the fused image. 

• Detect Changed Areas in the Fused Image Using the Improved FCM 

In reformulated fuzzy local information c-means clustering algorithm the local coefficient of variation is adopted 
to replace the spatial distance in FLICM which is defined as: 

C u =var(x)/(x) 2 (6) 

wherevar(x) and x are the intensity variance and the mean in a local window of the image, respectively. 

The value of C u reflects the gray-value homogeneity degree of the local window. 

The fuzzy factor for RFLICM can be defined as 

G'^N^+minacyCJ^Cu/a) 2 )] 

[X(l-u kj ) m llx r v k ll 2 ifC J u >= T u 

Or 

G' kl =S JENl [i/2-min((C J u /C u ) 2 ,(C u /C J l ,) 2 )] 
[X(l-u kj ) m l| Xj -v k ll 2 ifC J u <-C u (7) 

Where, v k is prototype value of k lh cluster, u k jis fuzzy membership of i th pixel with respect to cluster k, C\,is local 
coefficient of variation of neighbouring pixels, Cu is local coefficient of variation of central pixel and ~C u is the mean value 
of C J u that is located in a local window. 

The calculation of the membership partition matrix is performed as follows: 

U ^ =1/ { Zy =1 {[(ll X 1 - V kl| 2 + G 'ki)/(llx 1 -V J l| 2 + G ' J1 )] 1/(m " J) } } (8) 

Uk represents the fuzzy membership of the i lh pixel with respect to cluster k, and c is the number of clusters, 
llx r v k ll 2 is the Euclidean distance between object and the cluster centre. Here, the reformulated factor balances the 
membership value of the central pixel taking into account the local coefficient of variation, as well as the gray level of the 
neighbouring pixels. If there is a distinct difference between the results of the local coefficient of variation that are 
obtained by the neighbouring pixel and the central pixel, the weightings added of the neighbouring pixel in will be 
increased to suppress the influence of outlier. 

CONCLUSIONS 

In proposed method, a wavelet based image fusion scheme will be used to generate the difference image. 
The information of changed regions reflected by the mean-ratio image is relative in accordance with the real changed 
trends in multi temporal SAR images. On the other hand, the information of background obtained by the log-ratio image is 
relatively flat on account of the logarithmic transformation. Hence, complementary information from the mean-ratio image 
and the log-ratio image can be utilized to fuse a new difference image. Here, the RFLICM algorithm that incorporates both 
local spatial and gray information is proposed, which is relatively insensitive to probability statistics model. The RFLICM 
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algorithm introduces the reformulated factor as a local similarity measure to make a trade-off between image detail and 
noise. So, the proposed wavelet fusion strategy can integrate the advantages of the log-ratio operator and the mean-ratio 
operator and gain a better performance. The change detection results obtained by the RFLICM exhibited less spots than its 
pre-existence (i.e., FLICM) since it is able to incorporate the local information more exactly. 
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